Dynamics of cooperativity in chemical sensing among cell-surface receptors 
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Cooperative interactions among sensory receptors provide a general mechanism to increase the 
sensitivity of signal transduction. In particular, bacterial chemotaxis receptors interact cooperatively 
to produce an ultrasensitive response to chemoeffector concentrations. However, cooperativity be- 
tween receptors in large macromolecular complexes is necessarily based on local interactions and 
consequently is fundamentally connected to slowing of receptor conformational dynamics, which 
increases intrinsic noise. Therefore, it is not clear whether or under what conditions cooperativ- 
ity actually increases the precision of the concentration measurement. We explictly calculate the 
signal-to-noise ratio (SNR) for sensing a concentration change using a simple, Ising-type model of 
receptor-receptor interactions, generalized via scaling arguments, and find that the optimal SNR is 
always achieved by independent receptors. 
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In biological networks, cooperative interactions among 
components can sharpen input-output relations, increas- 
ing gain and enabling switch-like responses. The best- 
known example is the cooperative binding/release of oxy- 
gen by hemoglobin, which enables efficient transport of 
oxygen between the lungs and tissue. In sensory sys- 
tems, a well-studied example of cooperativity is receptor- 
receptor coupling in Escherichia coli chemotaxis |l[ . For 
sensing systems, there is an obvious advantage of high 
gain to amplify weak signals Q, particularly when com- 
bined with an adaptation system to broaden the dynamic 
range [3|. This advantage of high gain from receptor- 
receptor interactions raises the question: why has recep- 
tor coupling not evolved in other chemical sensing sys- 
tems, e.g. quorum sensing and eukaryotic chemotaxis? 

If the function of a sensory network is to reliably de- 
tect weak signals, then processing of the signal is only 
half the story. The other half is the suppression of 
noise and for weak signals the signal-to-noise ratio (SNR) 
generally governs information transmission and sensory 
performance Cellular signal transduction must con- 
tend with both noisy inputs (extrinsic noise), as well as 
the noise generated internally by the signal transduc- 
tion system itself (intrinsic noise). While the presence 
of signaling noise is well appreciated, the connection be- 
tween cooperativity and noise has received less atten- 
tion [H, Q . As we will show, receptor cooperativity and 
intrinsic signaling noise are inextricably linked via the 
statistical mechanics of receptors. The same cooperative 
interactions that give rise to enhanced sensitivity nec- 
essarily both amplify fluctuations and slow the rate of 
receptor-conformational switching, limiting both the re- 
sponse time and the ability of the system to reduce intrin- 
sic noise by time-averaging. Due to these tradeoffs, it is 
unclear when or whether receptor cooperativity actually 
increases sensory performance. 

The prevailing view of bacterial chemorcccptor opera- 
tion is that signal amplification from receptor coopera- 
tivity enables the chemotaxis network to reliably detect 



shallow gradients with exquisite sensitivity. In this Let- 
ter, we challenge this viewpoint by evaluating the strat- 
egy of using receptor cooperativity to enhance weak sig- 
nal detection in light of the tradeoffs between gain and 
intrinsic noise. Specifically, we calculate the SNR for a 
simple physical model of receptor-receptor interactions, a 
dynamical Ising-type model in which receptors have two 
conformational states, active and inactive, and neighbor- 
ing receptors prefer to be in the same conformational 
state (a.k.a. conformational spread [2]). The input to 
the network is the (changing) external ligand concentra- 
tion and the output is the time-averaged total receptor 
activity. This two-state model is a good description of 
bacterial chemoreceptors, in which the active and inac- 
tive conformational states are named for their relative 
abilities to activate a downstream kinase Q and the out- 
put of the network is the CheY-P concentration, which is 
effectively the time-average of receptor activity over the 
phosphorylatcd lifetime of individual CheY-P molecules. 

Importantly, our model is the first to both (i) incorpo- 
rate the dynamics of receptor switching, which is key to 
understanding the reduction in noise by time-averaging 
and (ii) base cooperativity on local interactions. Unlike 
the dynamic MWC models used in previous work [6|, we 
invoke no unphysical action-at-a-distance that allows the 
entire receptor cluster to flip conformational states si- 
multaneously. In fact, a recent study has provided direct 
evidence for local interactions underlying the ultrasensi- 
tive behavior of the bacterial motor [8j. 

We consider a cell with a total of N two-state receptors, 
which arc divided into m independent ID chains each 
of length n, with nearest-neighbor Ising couplings inside 
each chain, given by the Hamiltonian 
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where <jj = ±1 represents the active/inactive receptor 
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states, denotes nearest neighbors, and 



/ = Ae + log 



1 + [L]/K off 
1 + [L}/K on 



(2) 



is the free-energy difference between active and inactive 
states, which is a function of the ligand concentration 
[L], ligand dissociation constants K oS and K on , and the 
"offset energy" Ae = E on — E a g_, which is the energy dif- 
ference between active and inactive states in the absence 
of ligand 
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When ligand binding/ unbinding is much faster than re- 
ceptor conformational switching, the stochastic dynamics 
of switching is governed by a kinetic Ising model, where 
the occupation probablity p(u\ , . . .a n ;t) obeys the mas- 
ter equation 



d_ 

dt 



p{a x ,...a n ;t) = - p((Ti, . . .a n ;t) (3) 

+ '^2k i (-<T i )p(a 1 , . . . ,-cr,, . ..a n ;t). 



(For the case of slow ligand dynamics, see [4|.) We as- 
sume the switching rates k obey detailed balance and are 
local, i. e. depend only on the conformational states of the 
receptor and its nearest neighbors. Although the form of 
this dependence has not been measured experimentally, 
a simple, physically reasonable choice is Glauber dynam- 
ics [12J , which models interactions of the receptor system 
with a heat bath. With this choice, 
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7^0,-1 + cr l+ i) (1 - I3ai), (4) 



where a sets the intrinsic switching rate (typically 10 3 — 
10 4 /s), 7 = tanh(2J), and /3 = -tanh(//2) [H. 

We let the model system respond to a small step in 
ligand concentration (and associated step change in free 
energy Af ~ A log [L] ) over a given period of measure- 
ment, which we call T avg . We define the average activity 
of a receptor to be the probability the receptor is in the 
active state, with the total activity A being the sum of 
the individual receptor activities. The output of the sys- 
tem is the time-averaged change in total activity 14] . 
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X (t-t')Af(t')dt' (5) 
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where the dynamical susceptibility x relates changes in 
the average cluster activity to time-dependent changes 
in the free-energy difference and Eq. [6] defines the re- 
sponse function i?(r avg ). The system is assumed to be 
pre-adapted to the ambient ligand concentration, such 
that the free-energy difference between active and inac- 
tive states is zero prior to stimulation (i.e. the system is 
adapted to the most sensitive region of the input-output 
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FIG. 1. (A) Schematic diagram of cooperatively interacting 
receptors. Receptors interconvert between active and inactive 
states and preferentially bind ligand in the inactive state. The 
interaction energy J between neighbors favors receptors in 
the same activity state. (B) The normalized static response 
R(oo)/n for closed chains (rings) of various lengths n (for 
open chains, see II). (C) Response R, (D) noise a , and (E) 
SNR as functions of coupling strength J for a closed chain of 
n = 10 receptors for various averaging times r avg . In (E), the 
SNR is normalized by the averaging time. 



relation); this assumption appears to be correct for bac- 
terial chemoreceptors [9(. 

The sensitivity to signal is quantified by the static re- 
sponse R(oo) = 4 J °° x{t)dt, which, when normalized by 
chain length, is the factor by which the DC (i.e. infinite 
time) response of a coupled receptor is amplified rela- 
tive to that of an uncoupled receptor. The normalized 
static response R(oo)/n increases as a function of cou- 
pling strength from 1 to n, as shown in Fig.[Tj3. At zero 
coupling, J — 0, each receptor behaves independently 
and there is no cooperative amplification of the signal. 
As the coupling is increased, domains of adjacent recep- 
tors begin to effectively switch conformations together 
and the amplification R(po)/n is determined by the size 
of these domains. 

In practice, many cells have a limited measurement 
period and this can decrease the response i?(r avg ), as 
shown in Fig. [Tp. For a ID ring, this dependence takes 
the simple form 



-Rring(T aV g) = -Rring(oo) 1 
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where 



l/(a(l - tanh(2J))) 



(8) 



is the response time, which increases exponentially with 
coupling strength due to the well-known phenomenon of 
critical slowing down [lij]. Notice that for any finite av- 
eraging time, the response eventually falls to zero with 
increasing coupling strength because the system slows 
down and cannot respond to the input in the time avail- 
able. Intuitively, for large coupling J, the receptors be- 
come frozen in an all-active or all-inactive state, with 
switching between these states too slow to mediate a 
timely response to a changing input. 

The intrinsic noise, i.e. the variance in time-averaged 
activity, increases monotonically with coupling strength, 
as shown in Fig. [Tp. For the ID ring, 

^ingfovg) = 4 me (0)T(r avg /r^), (9) 



where r c is given in Eq. [8] and 



T(x) = 2{x + e- x -l)/x 2 



(10) 



For a given J, the maximal value of the noise occurs 
for zero averaging time (the "snapshot" limit), and is 
proportional to the static response as required by the 
fluctuation-dissipation theorem. Averaging for longer 
times substantially reduces the noise (via the factor 
7"(r a vg/r" ns )), but time-averaging becomes less effective 
with increasing J due again to critical slowing down, as 
the slower dynamics increases the correlation time r c , the 
time required for fluctuations in activity to decay. 

The relative uncertainty in sensing the concentration of 
ligand is determined by the signal-to- noise ratio (SNR), 
as shown in Fig. [IJD, normalized per receptor, 
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(11) 

Surprisingly, the optimal value of the coupling strength 
J to maximize the SNR is zero for all averaging times. 
That is, on a per receptor basis, independent receptors 
always have a lower total SNR than cooperative teams 
of receptors. For short averaging times, independent re- 
ceptors have both a larger response and lower noise than 
teams. For longer averaging times, cooperative teams 
have a larger response, but independent receptors still 
achieve a higher SNR because their rapid switching leads 
to more independent samples of receptor activity and 
thus to a much lower time-averaged noise. 

While the results in Fig. [Tp-E are shown only for 
n = 10, scaling analysis indicates that increasing the 
chain length n will never make cooperative teams favor- 
able with respect to uncoupled receptors. In the follow- 
ing, we derive how the SNR scales with the only two 
length scales in our model: the chain length n and the 
correlation length £, which is the length scale over which 



the conformational states of neighboring receptors are 
correlated in an infinite system. 

To best realize the cooperativity of a ID chain of 
length n, the coupling strength has to be set to J*, the 
smallest J that gives approximately maximal response 
(i?(oo) ps rt 2 ), and the averaging time has to be just 
long enough for R{T avg ) to approach this maximum, i.e. 
T avg ps t c (J*). Longer times will reduce the noise, but 
by exactly the same factor for both the cooperative team 
and independent receptors, oc l/r avg . 

For the cooperative team, r avg = r c ( J*) is long enough 
that R{r avg ) ps R(oo), while er 2 (T avg ) ~ c 2 (0) («-e. 
still roughly in the snapshot limit for noise). By the 
fluctuation-dissipation theorem, cr 2 (0) = i?(oo)/4, so 
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i? 2 (oo) R(oo) 



16ncr 2 (0) 
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*/4«£(.£)/4. (12) 



Thus the SNR for the cooperative team is approximately 
equal to the chain length. The last approximate equal- 
ity holds because the J* required for nearly maximal re- 
sponse yields an (infinite-chain) correlation length £(J*) 
comparable to the actual chain length n. 

By comparison, for n independent receptors, the same 
averaging time r avg = r c (J*) is long enough to reach 
the static response i?(r avg ) ps R(oo) = n, and for sub- 
stantial noise reduction by time-averaging, a 2 (r avg ) ps 
2CT 2 (0)/(ar avg ), so that 



SNR( J = 0) 
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Therefore, for cooperative teams to be favorable 
with respect to uncoupled receptors requires £(J^) > 
^ar c (J*), according to Eqs. IT21 and H"3l However, this 
will never happen for ID chains because the correlation 
time grows faster than the correlation length with in- 
creasing J in ID, specifically r c ~ £ 2 . In Fig. [2] we plot 
the ratio of SNRs for independent receptors versus co- 
operative receptors for closed chains. The ratio never 
favors cooperative receptors and, in the long-averaging- 
time limit, follows the expected scaling ~ £ for £ <C n 
and ~ £ 2 for £ 3> n. A shorter averaging time only 
makes things worse for cooperative receptors, as in this 
case the expected scaling is still ~ £ for £ <§; n, but ~ £ 4 
for £ > n [1]. 

Our results naturally generalize to higher-dimensional 
coupling. For 2D, the most natural topology for inter- 
acting membrane receptors, the ratio of uncoupled to 
coupled SNRs for long averaging times r avg ^> t c scales 
as SNR(J = 0)/SNR(J) -> nar c /R{oo) - « 
£°' 41 [1] [lij], as seen in Fig. [2] in the expected regime 
£ -C y/n. As in ID, independent receptors do better 
than coupled receptors, but this advantage grows more 
slowly with correlation length in 2D. 
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FIG. 2. Ratio of the SNR for independent receptors to the 
SNR for cooperative receptors as a function of the correlation 
length £ for a ID ring and a 2D square lattice of n = 100 
receptors with periodic boundary conditions, as well as mean 
field theory (MFT). Averaging times are ar avg = 100 (solid 
curve) and ar avg = oo (dashed curve), with power law plots 
for comparison (dotted lines). 



pendent receptors achieve SNR > 1 for the fewest recep- 
tors, and (ii) for a given number of receptors, indepen- 
dent receptors achieve SNR > 1 for the smallest stimulus 
strength. 

Our surprising result offers a fresh perspective on bac- 
terial chemotaxis, by indicating that the network is not 
simply optimizing SNR, since the observed receptor co- 
operativity in this system reduces the SNR by a factor 



of ~ 50 (18J. More generally, our result reveals that the 



benefits of cooperativity for sensing are far from obvious, 
potentially explaining the absence of receptor coopera- 
tivity in many sensory networks. 
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High dimensional (D > 4), or global coupling of recep- 
tors, which could be mediated by the cell membrane or 
by rapidly diffusing effectors, is described by the mean 
field limit of the Ising model. In this limit, the dynam- 
ics of cluster activity reduces to that of an overdampcd 
harmonic oscillator, and the normalized static response 
and correlation time are both equal, R MFT (oo)/n — 
ar MFT _ — ^j)^ where £ is the number of nearest 
neighbors. Consequently, for long averaging times, the 
ratio of SNRs for independent versus coupled receptors 
is constant. Thus, at best, increasing the dimension of 
receptor coupling yields parity between cooperative and 
independent receptors. 

More generally, the universal behavior of Ising mod- 
els near the phase transition and a rigorous bound on 
Ising critical exponents vz > 7 [17j suggest that critical 
slowing down is unavoidable and independent receptors 
will always optimize the SNR. We have also studied the 
addition of other potential noise sources, including slow 
ligand dynamics and static variation of receptor offset 
energies, and find that independent receptors still yield 
the best SNR [|. 

In summary, we developed a physical description of 
cooperativity based on the principle that allosteric in- 
teractions between receptor proteins are inherently local. 
From our simple Ising-type model, which encompasses a 
broad class of models, we elucidated the relationship be- 
tween cooperativity and intrinsic noise. We found that 
the slowing down of receptor switching due to coopera- 
tive interactions strongly impairs the SNR and that con- 
sequently the SNR is always highest for zero receptor 
cooperativity, even though the absolute sensitivity is op- 
timized for nonzero cooperativity. Since our SNR is nor- 
malized by receptor number and stimulus strength, our 
results show that (i) for a given small stimulus, inde- 
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